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Abstract 

A multistream model for spinless electrons in a relativistic quantum plasma is introduced by 
means of a suitable fluid-like version of the Klein-Gordon-Maxwell system. The one and two- 
stream cases are treated in detail. A new linear instability condition for two-stream quantum 
plasmas is obtained, generalizing the previously known non-relativistic results. In both the one 
and two-stream cases, steady-state solutions reduce the model to a set of coupled nonlinear ordinary 
differential equations, which can be numerically solved, yielding a manifold of nonlinear periodic 
and soliton structures. The validity conditions for the applicability of the model are addressed. 

PACS numbers: 52.27.Ny, 52.35.Qz, 52.35.Sb 
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I. INTRODUCTION 



The interest in relativistic quantum plasma systems is growing exponentially, not only 
because of the relevance to astrophysical problems but also due to the fast advances in 
strong laser-solid plasma interaction experiments. Indeed, the development of multi-Peta- 
Watt lasers will soon make possible to address simultaneous quantum and relativistic effects 
in laboratory plasmas Relativistic quantum kinetic models have been proposed, in the 
treatment of the plasma dispersion function for a Fermi-Dirac equilibrium [2|], and in the 
study of relativistic effects for quantum ion-acoustic wave propagation . Also, a covariant 
Wigner function theory for relativistic quantum plasmas described by the Dirac-Maxwell 
system has been suggested \^\, as well as the spinless (Klein-Gordon-Maxwell) analog has 
been presented jsl- The Klein-Gordon-Maxwell system of equations has been applied to the 
analysis of parametric scattering instabilities in relativistic laser-quantum plasma interac- 
tions [sl, while recent models have been introduced based on the Dirac-Maxwell equations 
describing the nonlinear propagation of light in Dirac matter j? | . 

Relativistic two-stream instabilities are traditionally known to be important for the elec- 
tron heating in intense laser-plasma interaction experiments js], as well as in astrophysical 



relativistic shocks 



9|], and could be important for pulsar glitches [lO|, 111, where superfluid 
neutrons and superconducting protons co-exist with relativistic electrons In addition, 
the two-stream instabilities between electrons and/or holes are also believed to exist in 
semiconductor plasmas |l3l-ll5|. 

In a lowest level of approximation than kinetic theory, quantum plasma hydrodynamic 
models are popular tools (see, e.g. 



16Nl9|). since they allow an efficient treatment of 



nonlinear phenomena, both from the analytical and numerical viewpoints. Starting from the 
Dirac-Maxwell system, there are hydrodynamic models for relativistic quantum plasmas 20| , 
which have been extended to incorporate particle-antiparticle effects in the wave propagation 
2l| . as well as relativistic quantum corrections to laser wakefield acceleration 22]. In these 
fluid formulations, we note, in particular, the non-trivial form of the relativistic extension 
of the quantum tunneling force in the momentum transport equation. Electromagnetic 
quantum hydrodynamic wave equations have also been considered including relativistically 
degenerate electron fluids 23|, based on previous spin-one-half hydrodynamic models 



24 



27l |. We note that hydrodynamic versions of the Klein-Gordon-Maxwell system have been 
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used in the past [28| , and recently in the context of laser physics [291] . 

In this paper, we introduce a relativistic multistream quantum plasma model starting 
from the Klein-Gordon-Maxwell^ystem of equations. We adapt the formalism according to 
the classical |3q| and quantum [31] multistream model for plasmas, and show its usefulness 
in a paradigmatic plasma problem, namely the linear and nonlinear features of the quantum 
two-stream instability in the relativistic regime. The simplicity of the Klein-Gordon equa- 
tions, in comparison to the Dirac equation, makes it a natural candidate for the extension of 
non-relativistic quantum plasma theories to the relativistic regime, when electron-one-half 
spin effects can be neglected. Therefore, a direct comparison to existing results on quantum 
plasmas can be obtained in an easier way. Moreover, an hydrodynamic formulation based 
on the Klein-Gordon-Maxwell system of equations strongly favors the development of new 
analytical and numerical tools for relativistic quantum plasmas. On the other hand, the 
domain of applicability is restricted to plasmas where the electron-one-half spin effects are 
not decisive, like in close to isotropic equilibrium configurations. For the parameters in 
this work, we have no quantized electromagnetic fields, so that quantum field theoretic re- 
sults involving e.g. pair creation are not included. Streaming instabilities in non-relativistic 
quantum plasmas are attracting considerable interest, since they display many surprising 
characteristics of pure quantum origin. Among these, we have a new instability branch for 
large wavenumbers, as well as new nonlinear spatially periodic solutions in the steady state 



case 321]. Besides formulating the relativistic version of the quantum multistream model for 



plasmas, the purpose of the present work is to extend the analysis of the quantum two-stream 
instability to the relativistic case. 

The manuscript is organized as follows. In Section II we introduce the multistream Klein- 
Gordon-Maxwell system of equations, casting it into a suitable fluid-like formulation. The 
one-stream case is treated in detail in Section III, where linear wave propagation is studied 
considering small amplitude perturbations around homogeneous equilibria, and a rich vari- 
ety of nonlinear periodic as well as soliton structures are found numerically. An existence 
criterion for solitary wave solutions is obtained. In Section IV, the relativistic quantum 
two-stream case is studied in depth. The linear relativistic quantum two-stream instability 
problem is fully characterized. In this regard, a main result of this work is the derivation 
of the instability condition in Eq. ([52D . which provides a natural generalization to the 
non-relativistic instability criterion |3l|. Both nonlinear periodic and soliton structures are 
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found numerically, and an existence criterion for localized solutions is obtained theoretically. 
Section V is dedicated to the final conclusions, including a detailed account on the validity 
domain of our model equations. 



II. THE KLEIN-GORDON-MAXWELL MULTISTREAM MODEL 

We here consider a relativistic multi-stream quantum plasma where the electrons are 
described by a statistical mixture of N pure states, with each wavefunction ipj satisfying the 
Klein-Gordon equation 

W^^j - c^V^^j - m^c% = 0, j = l,...,N, (1) 
where we have defined the energy and momentum operators respectively as 

W^th^ + e<l>, (2) 

and 

r = -ihV + eA . (3) 

Here, m and — e are the electron mass and charge, respectively, H is the Planck constant 
divided by 27i, c is the speed of light in vacuum, and 4> and A are the scalar and vector 
potentials, respectively. 

The electric charge and current densities are, respectively, 

N 



and 



p = E [rjWi^j + V'.ovv,)*] , (4) 



e ^ r 1 
The charge and current densities fulfill the continuity equation 

| + V.J = 0. (6) 

The self-consistent scalar and vector potentials are obtained from the inhomogeneous 
Maxwell's equations, using the Coulomb gauge V • A = 0, as 
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V20= -l(p + noe), (7) 

□A = M - ^V^ , (8) 

where a fixed neutrahzing ion background of the charge density eno was added, and where 
Eq and fio denote the vacuum electric permittivity and magnetic permeability, respectively. 
Here, the d'Alembert operator is 

1 (9^ 

The resulting Klein-Gordon-Maxwell system of equations ([T]) and ([7])-([HD describes the 
nonlinear interactions in relativistic quantum plasmas where spin effects and pair creation 
phenomena are negligible. In explicit form, the Klein-Gordon equation reads 



Following Takabayasi 



281], it is convenient to introduce a fluid-like formulation in terms 



of the eikonal decomposition 

ipj = Rj exp(iSj/h), (11) 

where the amplitude Rj and phase Sj are real functions. Separating the real and imaginary 
parts of the Klein-Gordon equation, we have 



1 / dSj \ /vTO A \2 2 2 fi^ORj 



C2 

and 



^.j (V5, + eA)^-mV = -^, (12) 



e 



2 dRj fdSj 

In terms of Rj and Sj, the charge and current densities are, respectively. 



mc^ ~^ \ dt 



and 
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N 



J = --E^?(V^.+eA). (15) 

Alternative hydrodynamic-like methods are also available, as for instance the Feshbach- 
Villars formalism 33|] where initially the Klein-Gordon equation is split into a pair of first- 
order in time partial differential equations. However, the resulting set of equations turns 
out to appear much more involved than in the present Takabayasi approach, which we use 
due to its formal simplicity. 

From now on we concentrate on the electrostatic case, and assume A = 0. In this 
situation, the relevant equations are 



1 fdS, 



dt 



- ec 



2„2 



2 dRj (dSj 



and 



- 2VRj ■ VSj = 0, 



(16) 
(17) 



e 

^^0 



2^ Qt 



no 



(18) 



The assumption A = is valid as long as the electric current is curl free, V x J = 0. This 
is true for longitudinal waves in a stationary plasma and parallel to a plasma beam. 



III. ONE-STREAM CASE 

We proceed next to study linear and nonlinear waves for the one-stream case (A^ = 1). 
For this case the sum in Eq. (ITS]) collapses to one term involving Ri = R and Si = S, and 
Eqs. (HSD-dlHD become 



and 



l_/dS 



- ec 



2„2 



h^UR 
R 



ecj) 



2VR-VS = 0, 



(19) 
(20) 



e 
£^0 



' R^ fdS 
mc^ \ dt 



-e(j)] + no 



(21) 
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A. Linear waves 



The system of equations (fT^ - (l?Ii) has the equihbrium solution 



R 



7 ' 



where 



7 



1 



pr, 



1/2 



(22) 



(23) 



is the relativistic 7 factor for a beam momentum p. Linearizing and assuming perturbations 
~ exp(z[K ■ r — Vtt]), where, for simplicity, we take K || p, and obtain 



(fi - Kvf = ^ + 



+ 



(24) 



ryo 4^2j^2 \ ^2 / 4^3^2^2 

where v = p/{'ym) is the equilibrium beam speed and Up = {noe^ /{meo)Y^'^. 

In the non-streaming limit J9 — > (and 7 = 1), Eq. (1241) is identical to Eq. (4.21) of 



Kowalenko et al. 



34| 



For the general one-stream case with waves propagating obliquely to the beam direction, 
the assumption A = fails, and one has to involve the full set of Maxwell's equations. 
However, in the one-stream case, one can start with the 3D dispersion relation in the beam 
frame, and then Lorentz transform the result to the laboratory frame. In the beam frame, 
the dispersion relation is 



l\2 



(A" 



i\2 



l\2 



+ 



i\2 



{K'Y - 



(25) 



where the primed Q' and K' are the angular frequency and wave vector of the plasma 
oscillations in the beam frame. 

To go from the beam frame to the laboratory frame, we assume for simplicity that the 
beam velocity is along the z axis. Then, the time and space variables are Lorentz transformed 
as t' = ■y{t — vz/c^), x' = x, y' = y and z' = 7(2; — vt). The corresponding frequency 
and wavenumber transformations are Q' = — vKz), K'^ = K^, K'y = Ky, and K'^ = 

^{K-^-vn/c^). 



The plasma frequency is transformed as uj' 



7. One easily verifies that the ex- 



pression {K') — ip!) /c = K — VL /c is Lorentz invariant. This yields immediately the 
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general dispersion relation for beam oscillations in the laboratory frame, 



(n - vK^f = ^ + 



+ 



(26) 



In the formal classical limit (/i = 0), we have from Eq. (12^ the Doppler shifted relativistic 



plasma oscillations 
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right-hand side of Eq. 



36j VL = Kv + u;p7 Using the limit Vt ^ Kv and i7 ^ in the 
one obtains 



{n - Kvf 



3 A'y^TTjP' 



(27) 



which is similar to the expression used by Serbeto et al. 37|] in the context of quantum free 



electron lasers, and where the last term in the right-hand side can be considered a quantum 
correction to the relativistic beam-plasma mode. 

On the other hand, in the non-relativistic limit c — > oo we have 7 = 1 and the familiar 
result [31 1 



{n - Kvf 



Am? ' 

describing Doppler-shifted quantum Langmuir waves. 
Normalizing into dimensionless units according to 



(28) 



(1' = ^ 



cK 



P 



one obtains (omitting the asterisks) 



mc 



H 



mc 



p_ 

,2 ' 



(29) 



where now 7 = {l+p^Y/"^ 
modes Vt = Q±, with 



l-t;2)-i/2. UHy^O and 



n') , (30) 
0, we solve Eq. fl5U]l to find the 



1 . 2 



m ^ V4 



1/2 



(31) 



It can be verified that both modes are stable > 0). It is reasonable to expand the last 



38 



39|, the 



result assuming a small H. Even for laser-compressed matter in the laboratory 
value of H presently does not significantly exceed 10~^. On the other hand, for conditions 



in the interior of white dwarf stars with a quantum coupling parameter exceeding unity we 
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reach the pair creation regime which can be safely treated only within the quantum 
field theory. The result is 



nl = ^^+2K'+ ^^{l - K') + 0{H') , (32) 
^l = l_^^^^{l_K'') + 0{H^). (33) 

Here fi+ is the pair branch which goes to infinity as — > 0. 

In dimensional units, the pair branch has a cutoff aX Vt = 2m(? /h. The is actually a 
backward wave (negative group velocity) for small wavenumbers, as mentioned by Kowalenko 
et al. y below their Eq. (4.23). 



B. Nonlinear stationary solutions 

Next, we consider nonlinear stationary solutions of Eqs. (IT^ - (PT]) in one spatial dimen- 
sion, of the form 

R = R{x) , S = --fmc^t + So{x), = 0(x) , (34) 

so that the original partial differential equation system is converted into a system of ordinary 
differential equations 



h'c'R" 



2 2 

p C 



2e'jmc 



,2/2 



R. 



RSq + 2R' Sq = 



-{jR'-no) + ^R\ 



(35) 
(36) 
(37) 



where the primes denote derivatives with respect to x. 

Equation fl36l) can be immediately integrated as R^S'q =constant. This relation is also 
equivalent to current continuity Jx =constant, which follows from the continuity equation (6) 



with dp/dt = 0. Assuming that R = ^Juq/'j and Sq = p where the plasma is at equilibrium, 
we have 



R^Sq 



noP 
7 



Sn 



nop 
7i?2 



(38) 
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which inserted into Eq. (155]) yields 



h^c'^R" + (fc^ + 2e7mc20 + e^"") R = . (39) 

Equations (|37ll and (1391) form a coupled nonlinear system for and i?, describing steady 
state solutions of our relativistic quantum plasma. 

Other special solutions (traveling wave solutions, alternative boundary conditions) could 
also be investigated, but we keep the above scheme, since then we can directly compare to 
the previous linear wave analysis. Indeed, Eqs. (1571) and (15^ admit the equilibrium 

i?2 = !!^, = 0, (40) 

7 

in the same way as the original Klein- Gordon- Maxwell system of equations (which also needs 
the equilibrium phase 5*). 

To proceed, we first transform into dimensionless variables according to 



R = ^= , 0=7 + — 2 , X = p = — , 5*0 = , 41 

^juq mc^ c mc mc^ 



so that the system for stationary waves becomes (omitting the asterisks) 



H'^R!' + (0^ - l)R 



V 

R^' 
R^6 - 1 



(42) 
(43) 



and 



^0 = ^, (44) 

where v = p/'y. 

In the system ( H2l) and ( H3l) . the variable x takes the role of a time-like variable, so that 
standard methods for ordinary differential equations can be applied. In this context, it is 
interesting to investigate the system around the equilibrium point if it admits only oscilla- 
tory (stable) solutions, or if it also admits exponentially growing (unstable) and decaying 
solutions. In the latter case, there is a possibility of finding localized solitary waves solutions 
with exponentially decaying flanks, which is not possible for stable cases. 
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Linearizing the system around the equihbrium (HUi) and supposing perturbations oc 
exp{iKx), we obtain the characteristic equation for the eigenvalues K 



1 = 0, 



which is the same as Eq. f l24l) with = 0. In the formal classical limit {H 
only linearly stable oscillations with 



(45) 

0), we have 



2„,2 



1 



(46) 



Hence, in the classical well-known, we do not have localized stationary solutions. 

In the quantum case {H ^ 0) the situation is more complex. The characteristic equation 
can be solved yielding K'^ = K'^, with 



7P 



— + 1P 



nl/2> 



(47) 



Expanding for small H, we obtain 



(48) 
(49) 



Note that is the quantum extension of the branch in Eq. fHBl) . while has no classical 
analog. 

To investigate the stability of the equilibrium, it is useful to rewrite the characteristic 
equation as 



7 



(50) 



which gives a second degree equation for K^. Since the characteristic function F{K'^) satisfies 



7 



4p2 



^(0) = r^>i, F\K'>^\<^ 



p^ 



(51) 



and has a pole at K"^ = Ap^ / H'^ 



lim FiK'^) = ±oo , 



(52) 
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we will have stable oscillations provided that the minimum value F^-m in the branch < 
K"^ < Ap^/H'^ satisfies 

i^min < 1 . (53) 

In this case F^K"^) intercept the value 1 at two positive values, so that the characteristic 
equation has only real solutions. The situation is summarized in Fig. 1. 



F F 




FIG. 1: on the left, a typical case of stable linear one-stream oscillations when -Fmin < 1- On the 
right, a typical unstable case where -Fmin > 1- The pole is at = Ap'^/H'^ 

Working out Eq. (1531) we find that 

Further analysis shows that the condition (jSlD can be written as 

H < //max = 27(7'/' - 7"'/') , (55) 

as the final condition for stable linear oscillations. It follows that the necessary existence 
criterion for nonlinear localized (soliton) stationary solutions is H > i/max- In Fig- 121 a 
graph of -ffmax is plotted as a function of the velocity v (measured in units of c), which 
shows that the quantum range for stable oscillations is increased for increasing relativistic 
effects. In the non-relativist ic limit p ^ 1, we have H < = a,s the condition for stable 
oscillations, or in dimensional units, hwp < mv"^. 

For the nonlinear system fH2]) and (1431) . a Hamiltonian form can be obtained with the 
further transformation 

R-i- iR, X ^ X , (56) 

H 

so that 
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10 



8 - 

6 - / - 

^max / 

4 - 

2 y 

0.2 0.4 0.6 0.8 1 

V 

FIG. 2: maximum quantum parameter i?max in Eq. (j55p as a function of the velocity v measured 
in units of c 



where 



Since we arrive at an autonomous Hamiltonian system, one has the energy integral 

+ ^ + (59) 



2 2 

Restoring dimensional variables, we have the conserved quantity 



"-/-27, (60) 



or 



mc^ dx 



+ 



2e0 R 



— 7 + 
^0 V 



1 /rfi?' 



i?2 ' 



(61) 



where /3 = f/c. The obtained conservation law can be used to verify the accuracy of 
numerical simulations. 
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FIG. 3: Spatial profiles of (p, R, and the electron number density Ng = B^cf) (top to bottom panels) 
for 7 = 1.3, and H = 0.01 (left column), H = 0.1 (middle column) and H = 0.2 (right column). 
The solution was set to (^(0) = 1.15, ii(0) = 0.98 and (^'(0) = R'{Q) = at the left boundary 



Numerical solutions of the nonlinear system fH2]) and (1431) are presented in Figs. |3HS1 For 
the non-localized solutions in Fig. [31 initial values on 0, i?, and their first derivatives were set 
on the left boundary and solution was integrated using the standard 4th-order Runge-Kutta 
method. For the localized solutions in Figs. HHHl boundary conditions on cj) and R were 
fixed on both the left and right boundaries, and the solutions were found with iterations 
based on Newton's method. Figure [3] shows large amplitude oscillations for 7 = 1.3 and 
different values of if, such that small-amplitude oscillations are linearly stable, in the sense 
discussed above. We see that there is one short and one long length-scale, corresponding to 
Kj^ and for the linear oscillations in Eq. ( H71) . Here the small-scale oscillations are due 
to the quantum diffraction effect, while the large-scale oscillations are related to wakefield 



oscillations which are well-known in classical plasmas [41i- When these two length-scale 
become comparable, i.e. for large if, the two length-scales interact and the oscillations 
become more irregular. 

For parameters where the oscillations are exponentially decaying or increasing (unstable 
oscillations), we have the possibility of localized solutions in the form of dark or grey solitons. 
It turns out that the coupled system of equations ( l42l) and ( l43i) supports a wide variety of 
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FIG. 4: Spatial profiles of (j), R, and the electron density N^, = R'^cj) (top to bottom panels), showing 
single, double and triple dark solitary waves (left to right columns) for the zero beam speed case 
V = with H = 0.01. The solution is set to (p = \R\ = 1 at the left and right boundaries {R = — 1 
on the left boundary for the single and triple dark solitons) 

nonlinear localized structures. Due to quantum diffraction effects, the plasma can develop 
dark solitary waves with one or more electron density minima. In Fig. HJ we see different 
classes of dark solitary waves for the case when the plasma is at rest, v = 0, and H = 0.01. 
Since p = v = 0,m this case, the term proportional to in the right-hand side of Eq. (142 p 
vanishes, and R can continuously go between positive and negative values. We see single, 
double and triple dark solitons, where R is shifted 180 degrees (from negative to positive) 
on between the left and right sides of the single and triple dark solitons. The dark soliton 
with single electron minimum is the same type as found in Ref. 421] for a non-relativistic 
quantum plasma. The solutions with multiple density minima are somewhat similar in shape 
to the multiple-hum p optic al solitons predicted in relativistic laser-plasma interactions in 



the classical regime 



43 



441 . On the other hand, in Figs. and El we consider solitons in a 



streaming plasma with finite speed f > 0. We recall that the existence condition for solitons 



is H > 27(7 



1/2. 



■7 



[where 7 = 1/(1 — f^)^/^], which puts an upper limit on v for a given 



value of H. For example, for H = 0.01, shown in Fig. [5l we have w < 0.1, while for H = 0.5, 
shown in Fig. El we have v < 0.58 for solitons to exist. For H <^ 1 and non-relativistic 
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FIG. 5: spatial profiles of (j), R, = R^cj)-, and (5o —px)/H (top to bottom panels), for H = 0.01, 
and V = 0.01 (left column), v = 0.05 (middle column), and v = 0.09 (right column). The solution 
was set to (/) = 7 and R = 1/^/7 at the left and right boundaries. We see grey solitons with 
non-zero electron density in the center. Bottom panels show the phase shift (5o — px)/H of the 
wavefunction 

f <^ 1, the existence condition for solitary structures becomes < H, or in dimensional 
units, mt>^ < fiUp. A general feature of the propagating solitons is that the electron density 
is non-zero at the center of the soliton, hence they are grey solitons. Furthermore, as the 
speed increases, the amplitudes of the solitons decrease and their tails become oscillatory 
when the speed approaches the maximum allowed speed, as can be seen in the right-hand 
columns of Figs. O and O There is also a complex phase shift proportional to 5*0 in the 
total wave function ip due to the relation (155]) . The plot of {Sq — px)/H in Fig. shows 
how the phase (in radians) is shifted between the two sides of the solitons. As f — 0, the 
phase jumps abruptly a value of ~ vr in the center of the soliton, while solitons with higher 
speeds have smaller and smoother jumps in the phase, as can be seen in the bottom panels 
of Figs. and El 



16 



1.5 



-10 
2 



CC 1 



0^ 
-10 



-10 
X 4 
^ 2 

Q- 
I 

- 



03 



-10 



v=0.1 



1.5 



10 -10 
2 



10 -°10 



10 -10 
4 



2 


10 -10 



v=0.3 



1.4 
1.2 
1 



10 -10 
2 



10 -10 
2 

— 1 



10 -10 
4 

— 2 


10 -10 



v=0.55 



10 



10 



10 



10 



FIG. 6: Spatial profiles of cp, R, = R'^cp, and {Sq — px)/H (top to bottom panels), for H = 0.5, 
and V = 0.1 (left column), v = 0.3 (middle column), and v = 0.55 (right column). The solution was 
set to (/) = 7 and R = 1/^/7 at the left and right boundaries. We see grey solitons with non-zero 
electron density in the center. Bottom panels show the phase shift (Sq —px)/H of the wavefunction 

IV. THE TWO-STREAM CASE 

We next consider linear and nonlinear waves for the two-stream case (A^ = 2). For this 
case, stream is represented by a wavefunction ipj = Rj exp{iSj/h), j = 1, 2, and we have the 
Klein-Gordon-Poisson system of equations 



^ ^^^^ eA' -{VS,r-m'c' = ^^ 
J 



c2 V dt 

a, (us, - 4|) + ifi m 

\ at at \ ot 



2VRj ■ VSj = 



e 

£0 fnc"^ 



'dSj 
'dt 



e(f)\ +no 



(62) 
(63) 

(64) 



which describe a relativistic quantum two-stream plasma in the electrostatic approximation, 
using physical variables. 
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A. Linear waves 



Similar to the one-stream case, we have the equihbrium 



Si = —'jmc^t + p ■ r , 5*2 = —jmc^t — p • r , (65) 

for two symmetric counter-propagating electron streams. 

Linearizing the governing equations and assuming plane-wave perturbations with the 
wavenumber K || p and the angular frequency Q, we obtain the dispersion relation 

F{n) = 1 , (66) 
with the characteristic function F{Q) defined by 

F(n) = ^S- 4m^c^ - h\n^ - c^K^) 

where 

a;,= (^^j ' - = P/(7-). (68) 
In the classical limit, viz. = 0, we obtain the same results as in the description of 

n 

classical cold relativistic electron beams using the fluid theory [8|. It can be treated in full 
analytical detail. We have the dispersion relation 

Equivalently, it is useful to write the classical dispersion relation as Fci^) = 1, with the 
characteristic function 



(70) 



1 

+ 



obtained by setting ^ = in Eq. ( p7l) . The dispersion relation turns out to be a quadratic 
equation for Q"^, hence Fc{^) should attain the unity value four times to prevent instability. 
Graphically (see Fig. [7]) we conclude that 
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FciO) < 1 (71) 
is the condition for linear stability. This shows that the wave-numbers such that 



(72) 



are linearly stable. The same conclusion is reached analyzing the potentially unstable mode 
in Eq. (l69ll . In comparison with the non-relativistic stability condition > Up/v'^, we 
note that the relativistic effects are stabilizing, since they imply a smaller unstable range in 
wave-number space. 

Fc[£2] 




-Kv Kv 




FIG. 7: On the left, stable linear waves satisfying Fc{0) < 1, with the non-quantum characteristic 
function Fc{^) given by Eq. (f7D|) . On the right, unstable linear waves 



Setting fl = ifli for real fli and using Eq. 



we obtain 



273/2 



(73) 



as the maximum growth rate, which also becomes smaller due to relativistic effects. 

On the other hand, in the quantum but non-relativistic {h 7^ 0, and c 00) limit we 
have 



E 



0. 



(74) 



We will not discuss the non-relativistic case, since this has been already done in the 
past [31]. The non-relativistic well as the non-quantum case, can be solved in full 

analytical detail, because in both situations the dispersion relation is equivalent to a second 
degree polynomial equation for D?. 

We now turn our attention to the fully quantum- relativistic dispersion relation (1^ . Due 
to the symmetry, we can restrict the treatment to positive frequencies, wave-numbers and 
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beam velocities. Equation (1^ is equivalent to a fourth degree polynomial equation for Q'^, 
which can be analytically solved in terms of cumbersome expressions, or solved numerically. 
However, it is more informative to first analyze the behavior of the characteristic function 
F{Q) in Eq. (jSZD, which is mainly determined by the poles at 



h ~ 

,2 



+ 



'fmc 



'ymc 



h 



+ 



+ 



h 

h 
h ' 



Kv 



1/2 



r 



+ Kv\ + 



1/2 



—2 



Kv 



h 



7^ 



+ Kv\ + 



1/2 



1/2 



(75) 
(76) 
(77) 
(78) 



paying attention just to the positive values. More precisely, i^i > provided hK < 2p, 
otherwise the positive pole is at —Qi. It can be shown that one has the ordering 



< < ^3 < ^^4- (79) 

Here, Qi and ^2 have classical counterparts as /i — > 0, while and ^4 are associated with 
pair branches, without classical counterparts. 

We note that the case hK = 2p is degenerate, since then one has i^i = and the 
dispersion relation becomes a third degree polynomial equation for Q^. The solutions to this 
particular case always correspond to (marginally) stable modes, not considered any further 
here. 

A tedious analysis shows that the characteristic function has the following properties 



lim^ F{Q) = TOO , lim F{Q) = ±00 , 

hK <2p^ F{ni<n<n2) <o (80) 

SignF(O) = -SignF"(0) = -Sign{hK - 2p) . 

Moreover, F{Q) tend monotonously to zero as i7 — > 00. These results imply that the 
wave-numbers satisfying hK > 2p are always stable, since in this case the characteristic 
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function has the topology shown in Fig. [HI where F{Q) always intercepts the value unity 
four times. 

F[Q] 




Q 



FIG. 8: Generic behavior of the characteristic function F{i}) in Eq. (I67p for hK > 2p. Since F = 1 
at four positive frequencies, the corresponding wavenumber is stable 

On the other hand, the case hK < 2p is potentially unstable, according to the minimum 
value F{0). If F{0) < 1, the characteristic function attains the unity value at four positive 
frequencies, corresponding to four linearly stable waves. Otherwise, when -F(O) > 1 there is 
a (purely imaginary) solution for the dispersion relation, and hence instability. The whole 
scenario is summarized in Figs. M and [TDi 



F[Q] 



1 
















£22 




£23 




£24 



FIG. 9: Generic behavior of the characteristic function F{^}) in Eq. (|67|) for HK < 2p, in the stable 
cases where F{0) < 1 
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F[Q] 




FIG. 10: Generic behavior of the characteristic function F{Q) in Eq. (j67p for hK < 2p, in the 
unstable cases where F{0) > 1 



In summary, besides HK < 2p we have 



F(0) 



1 + h^Ky{4m\^) 



> 1, 



(81) 



as a necessary condition for unstable linear wave propagation in our two-stream relativistic 
quantum plasma. Rearranging the instability conditions found, we combine them according 
to 



P 



(82) 



It can be verified that Eq. (1521) reproduces the non-relativistic results |31 |. 

The instability condition (I82p . together with a numerical solution of Eq. are depicted 
in Fig. [m Here we assumed that = fi^ + zfi,, where fi,- is the real frequency and fij the 
growth rate, and plotted Qi/up as a function of = h^ujp/rn?v^ and i^2 = K'^v^ juJ^ for 
different values of 7. The case 7 = 1 corresponds to Fig. 1 of Ref. 31|, while 7 > 1 show the 
relativistic effects on the instability region. We note that in the formal classical limit the 
largest unstable wavenumber becomes smaller as 7 — > cxd. On the other hand, the height of 
the upper curve in the instability diagram scales as 72, so that in this sense the combined 
quantum-relativistic effects tend to enlarge the unstable area. Ultimately, however, quantum 
effects stabilize the sufficiently small wave numbers, no matter the strength of relativistic 
effects. An interesting quantum effect is the appearance of an instability region at large 
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FIG. 11: Unstable function of = }T?ujp/m'^v^ and = K'^v'^ /J^ for different values 

of 7, where the unstable cases, given by Eq. (|82p are between the lower and upper boundaries of the 
colored region. The color corresponds to the normalized growth rate 0,i/ujp obtained numerically 
from Eg. (I66|) . The case 7 = 1 (upper left panel) is the non-relativistic case corresponding to Fig. 1 



of Ref. 



3l[ |. For increasing values of 7, the instability region is shifted towards larger values of 



and the growth rate decreases 

wavenumbers = K^v'^/Up for moderately small values of = h'^Up/m'^v'^, which does 
not have a classical counterpart. Finally, it should be noted that simultaneously > 1 
and 7 > 1 in Fig. [TT] correspond to extremely high electron number densities, comparable 
to those in the interiors of white dwarf stars and similar astrophysical objects. 

B. Nonlinear stationary solutions 

We consider the one-dimensional version of the system (I5^ -( 1M|) and stationary solutions 
of the form 



Ri,2 = Ri,2{x) , 5*1,2 = -'-ymc^t + cri,2(a;) , 4> = 4>{x) 
Equations (!^ are then equivalent to 
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(83) 



^(^M) = ^(i?X) = 0, (84) 



where the primes denote x— derivatives. Assuming 



Rl = Rl = ^, a[ = -a',=p (85) 



at equihbrium, applying the transformation 



Ri,2 = ^= > 0=7 + — ^ , ^ = , (86) 



and using Eq. to ehminate a[ ^ from Eq. (152]) , we readily derive the system of equations 
(omitting the asterisks) 



l)i?i 



R'K!, + (0' - l)i?2 



4i?3 



4i?3 



(87) 
(88) 



and 



0" = (i?? + i?2)0_i^ (89) 

which predict nonlinear stationary solutions of a relativistic quantum two-stream plasma. 
Hereif and v are defined as in the one-stream case. 

Linearizing around -R12 = 1/(27) ^■iid = 7, and supposing perturbations oc exp(z_ft'x), 
we obtain a quadratic equation for ^ which is also obtained from Eq. ( |66l) setting = 0. 
Proceeding as before, we formally obtain the same existence condition as Eq. (1551) for 
periodic solutions. The two-stream nonlinear solutions can be constructed from the one- 
stream cases by using = |_R2| = |-R|/"\/2, where R is obtained by solving the system 
fH2]) - fH51) . The signs of R\ and Ri are arbitrary. In Fig. [T21 we show a numerical solution 
of the system (1621) -(l89l). where the profiles of and are identical to the ones in Fig. [3l 
and with Ri = R2 = Rj V^. In Fig. [13], we perturbed this solution by using different values 
of -R'i(O) and R'2{0) at the left boundary. The general behavior of the solution in Fig. [T5] 
remains similar as in Fig. [T21 but differences in the two solutions can be seen in the details. 



24 



H=0.01 



H=0.1 



0.5 



10 20 







0.5 



10 20 

1.5 



10 



1.5 



0.5 



0.5 

20 
2 
1.5 



I ^ ' 0.5 L 

3 10 20 



0.5 

10 20 °0 

4 



20 



10 



H=0.2 




10 20 




10 20 




10 20 



AlllililllBI 



10 20 



FIG. 12: Spatial profiles of (p, Ri, R2, and the electron density = {Rf + i?!)*/* (toP to bottom 
panels) for H = 0.01 (left column), H = 0.1 (middle column) and H = 0.2 (right column). The 
solution was set to 0(0) = 1.15, Ri{0) = i?2(0) = 0.98/^/2 and 0'(O) = R[{0) = i2^(0) = at the 
left boundary 

Stationary localized solutions are shown in Figs. [14] and [H] for the non-streaming {v = 0) 
and streaming {v > 0) cases, respectively. In both cases, we found only localized solutions 
corresponding to |-Ri| = I-R2I = \R\/V^ (where R is the one-stream solution), but no other, 
more complicated cases. For the non-streaming solutions in Fig. [T3] we show an example 
with |_Ri| = |_R2| = \R\/V^ in the first column (with opposite signs on Ri and i?2)- When 
the solution was forced to an anti-symmetric Ri and a symmetric R2 in space, the numerical 
solution converged to solutions where either Ri or R2 took the shape of a one-stream dark 
soliton, while the other part tended to zero (or as small as possible) close to the soliton. 
For the streaming case in Fig. [T51 we also only found localized solutions corresponding to 
\R^\ = \R^\ = \R\/y/2. Similar as in the one-stream case, we have a maximum beam speed 
of f ~ 0.1 for the existence of localized solutions, and as the beam speed approaches this 
value, the amplitude of the soliton decreases and becomes oscillatory in space. 

One further issue is the stability of the localized solutions in the streaming cases f > 0. 
Far away from the localized solution, the plasma can be considered to be homogeneous, and 
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FIG. 13: The same as in Fig. [H but using R[{0) = 10 and i?^(0) = -10 (left column), R[{0) = 1 
and R'2i0) = -1 (middle column), and R[{0) =0.1 and i?^(0) = -0.1 (right column) 



one can perturb the equilibrium and study plane wave solutions proportional to exp{—iQt + 
iKx) for real- valued K and complex-valued Q with unstable solutions if the imaginary part 
of Q is positive. In the one-stream case, studied in Section 2, all solutions were found 
to be stable in time, while in the two-stream case, studied in Section 4, we have K for 
which the solutions are unstable. Hence, for the two-stream case the system is sensitive 
to perturbations far away from the localized structure. The general stability analysis for 
localized solutions can be carried out with normal mode analysis by perturbing the nonlinear 
equilibrium solution of the system (1^ - (!M1) as Ri^2{x,t) = Ri^2{x) + Ri{x) exp{—iQt), 
Si,2{x,t) = —'jmc^t + ai^2{x) + Si^2{x) exp{—iQt) and 4>{x,t) = 0(x) -|- 0(x) exp{—iQt), and 
assuming that the perturbed quantities vanish at |x| = oo. This leads to a linear eigenvalue 
problem with eigenf unctions -Ri(x), Si^2{x) and 0(x), and eigenvalue fl. Solutions with fl 
having positive imaginary parts are unstable and will grow exponentially with time. 
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FIG. 14: Spatial profiles of 0, Ri, R2, and A'^e = + .Ri)*^ (top to bottom panels) for the zero 
beam speed case v = with H = 0.01. The solution is set to (j) = \Ri\ = |i?2| = l/^/2J at the left 
and right boundaries 

V. SUMMARY AND CONCLUSION 



In this paper, we have presented a multistream model for a relativistic quantum plasma, 
using the Klein-Gordon model for the electrons. We have treated the one- and two-stream 
cases in detail. We have derived dispersion relations for the linear beam-plasma interactions 
in the one-stream case, and for the streaming instability in the two-stream case. The system 
exhibits both plasma oscillations close to the plasma wave frequency that reduce to the 
Langmuir oscillation frequency in the classical limit h ^ 0, and pair branches which do not 
have a classical analog. Also, there is a new instability branch for large wavenumbers of 
pure quantum origin. A main result of this work is the derivation of the instability condi- 
tion in Eg. ( l82l) . which provides a natural generalization to the non-relativistic instability 
criterion [Sll]- Another important result is the condition (ISSl) for the existence of periodic 
(stable) oscillations and exponentially growing and decaying (unstable) steady-state oscilla- 
tions, where the latter exist only below a given electron beam speed. Similar to the classical 
plasma case 45|] , this furnishes an existence condition for periodic solutions, which exist only 
for cases with exponentially growing and decaying solutions. A rich variety of nonlinear 
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FIG. 15: Spatial profiles of i?2, and N^, = {R\ + R'2)(j) (top to bottom panels), for H = 0.01, 

and V = 0.01 (left column), v = 0.05 (middle column), and v = 0.09 (right column). The solution 
was set to = 7 and Ri = R2 = Xj^flF) at the left and right boundaries. We see grey solitons 
with non-zero electron density in the center 

solutions has been numerically found, including solitary waves with one or more electron 
density minima and an associated positive potential. It has been noted that the amplitude 
of the solitons decreased for increasing beam speeds, with increasingly oscillatory tails as 
the beam speed approached its maximum value for the existence of soliton solutions. 

Our model can be applied to situations where the quantum statistical thermal and elec- 
tron degeneracy pressure effects are small. The relative importance of these two effects can 
he characterized by the degeneracy parameter x = Tp/T, where Tp = fi?{?>Tx'^nQY^^ / {2KBm) 
is the Fermi electron temperature, and is the Boltzmann constant. When x < 1, the 
thermal pressure dominates, while when x > 1, the degeneracy pressure dominates. Our 
model is applicable when kbT <^ mv'^ for x > 1 and kbTf <^ m-v'^ for x < 1. Strong cou- 
pling (collisional) effects can be neglected when the coupling constant F = e^rig^^/ {A-ksoKbT) 

1 /3 

(for X < 1) or F = e'^riQ /{AttsoKbTp) (for x > 1) is small. For x > 1? the Pauli blocking 
further helps to reduce the effect of collisions [46,]. The pair creation phenomena have been 
neglected here, because our model excludes quantized fields. Hence, fkVp needs to be much 
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smaller than 2mc^ 40| . Working out the weak coupling and no quantized field assumptions, 
formulated as F < 1 and hup/ {2mc^) < 1, we have (using SI units) 



logioT>^logioTio-4.8, if x<l, (90) 



logiono>28.8, if x>l, (91) 

and 

logiorio < 38.9. (92) 

as the condition for the applicability of our theoretical model. 

Finally, it should be noted that our investigation neglects electron-one-half spin effects, 
which are justified since we used an unmagnetized quantum plasma model, and hence there 
are no spin couplings to a magnetic field. In conclusion, we stress that the present investi- 
gation of linear and nonlinear effects dealing with relativistic electron beams in a quantum 
plasma is relevant for high intensity laser-plasma interaction experiments white dwarf 



stars 471, |48|], and neutron stars 10l-ll2l|. where both quantum and relativistic effects could 



be important. 
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